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Abstract 



We study the dynamics of a nearly incompressible fluid via (classical) effective field 
theory. In the kinematical regime corresponding to near incompressibility (small fluid 
velocities and accelerations), compressional modes are, by definition, difficult to excite, 
and can be dealt with perturbatively. We systematically outline the corresponding 
perturbative expansion, which can be thought of as an expansion in the ratio of fluid 
velocity and speed of sound. Our analysis is fully relativistic, but, if desired, the 
non-relativistic limit can be taken at any stage. We demonstrate the potential of 
our method by some sample calculations: the relativistic corrections to the emission 
and to the scattering of sound by vorticose (possibly turbulent) motions, and, more 
interestingly, the vortex- vortex long range interaction mediated by sound. These results 
are, to our knowledge, new. For the last one, we find a 1/r^ potential between localized 
vortex configurations (e.g., vortex rings), whose strength and detailed structure depend 
on a certain multipole moment of the vortices under consideration. We also improve 
on the so-called vortex filament model, by providing a local field theory describing 
the dynamics of vortex-line systems and their interaction with sound, to all orders in 
perturbation theory. 

1 Introduction 

Hydrodynamics is notoriously hard to solve. This is usually blamed on the non-linear na- 
ture of its equations of motion, and we will not quarrel with that. It is instructive how- 
ever to ponder what makes hydrodynamics stand out with respect to other non-linear field 
theories (general relativity, for instance), where perturbative techniques exist that can be 
applied efficiently to a variety of physical situations. Our position is that there are certain 
configurations — the vortices — that in general cannot be dealt with in perturbation theory. 
This is discussed at length in ref. [T], where the quantization of hydrodynamics was at- 
tempted, but the problem is there already at the classical level. 

Consider, for instance, a non-relativistic incompressible fiuid, where the only allowed 
configurations are vortices, i.e., divergence-free velocity fields. In that limit, the continuity 
equation is just the constraint V ■ v = 0, the only degree of freedom is the vorticity field 
u = V X V, and the Euler equation can be rewritten as 
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If one neglects the non-linearity on the right hand side there is no time-evolution. The 
dynamics are completely dominated by non-linearities. 

An equivalent viewpoint is the effective field theory one pLj: Vortices correspond to ex- 
citations that, because of an infinite-dimensional symmetry (volume-preserving diffeomor- 
phisms), have vanishing energy in the limit of vanishing velocities, regardless of the excita- 
tions' spatial gradients. It is then easy for the system to excite these zero-energy gradients, 
and to make them exit the regime of validity of perturbation theory. Indeed, in field theory 
gradient energy counts as quadratic potential energy in the Hamiltonian. In the absence of 
a quadratic potential, the dynamics are dominated by anharmonicities. 

The situation is strikingly different for compressional modes (in a compressible fluid). 
At the equation of motion level they correspond to perturbations 5p in the density field, 
obeying a wave equation plus non-linear corrections: 



Even neglecting the non-linearities, there are some non-trivial dynamics corresponding to 
the free propagation of sound waves. Then — for small 5p — the non-linearities can be treated 
as small corrections to this. At the effective field theory level, compressional modes are 
standard gapless fields, which carry kinetic and gradient energy, and for which non-linear 
terms correspond to standard perturbative interactions p. 

These considerations indicate that there should exist physical situations in which com- 
pressional modes can be dealt with in perturbation theory, even though there might be an 
underlying vorticose 'background' that cannot. A nearly incompressible fluid provides such 
a system; there, by definition, compressional modes are difficult to excite. 

It should be remembered that near incompressibility is not an intrinsic property of certain 
fluids, but rather a kinematical regime that exists for all fluids: for fluid flows that are much 
slower than sound, any fluid behaves as nearly incompressible [2], and vice versa, for fluid 
flows that are as fast as sound, any fluid is quite compressible. The analysis that follows 
thus applies to any fluid in the appropriate kinematical regime. 

We will develop an effective field theory of sound interacting with vortices, and we will 
outline the systematics of the perturbative expansion that can be carried out close to the 
incompressible limit (sect. [2])j^ Some physical aspects of our analysis are not entirely new. 
For instance, the study of sound generation due to fluid flow (as opposed to time-dependent 
boundary conditions) was pioneered in the 1950's by Lighthill [21 H]- Since then, its study 
has had a long history full of many contributions and contributors, the details of which can 
be found in the relatively recent texts [5l [6] . 

The basic idea of Lighthill is to identify the non-linear terms in (|2]) (the 'dots') as the 
source of sound, and to try to solve that equation perturbatively. In a sense our field theory is 
doing just that. However, by insisting on field theoretical ideas formulated in terms of a local 
action we can be completely systematic in how we carry out such a perturbative expansion — 
simply because perturbation theory in quantum field theory has been exhaustively studied. 

^We will ignore dissipative effects due to viscosity and heat conduction throughout our paper. This is 
consistent at low frequencies and long wavelengths, because these effects correspond to higher order terms 
in the derivative expansion. An initial attempt at incorporating dissipation in our effective field-theory 
language can be found in . 
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We can organize each computation in terms of Feynman diagrams. By restricting to tree-level 
diagrams, we will be, in fact, just solving the classical equations of motion perturbatively. 
The Feynman diagram language is an extremely powerful organizational and computational 
tool to do that. 

It should be noted that Lighthill's approach to eq. ^ is the hydrodynamical analog of 
the so-called post-Netwonian expansion of general relativity (see e.g. [7j), which has been 
recently recast into an effective field theory language by Goldberger and Rothstein P]. Our 
paper owes much to theirs. 

Our first application of our tools will be to reproduce two preexisting results in the 
literature: the rate at which vorticose motions emit sound (sect.jsj), and the cross-section they 
have for scattering it (sect. |4]). However, since our field theory is relativistic by construction, 
we will be able to provide the relativistic corrections to these results at no additional cost. 
To the best of our knowledge, these have never been computed before. 

On top of providing powerful computational tools, the effective field theory also offers a 
novel viewpoint on the dynamics, which can lead to simple predictions of previously over- 
looked phenomena. For instance, since vortices interact with sound waves, and since these 
are gapless, from the field theoretical viewpoint it is obvious that sound waves can mediate 
long-range interactions between vortices. These are in addition to the well known purely 
kinematical "dragging" interactions, whereby the long-distance tail in each vortex's velocity 
profile drags all the other vortices with it. We are not aware of any mention of our sound- 
mediated interactions in the literature. We will compute the leading contribution to the 
associated potential (sect.jsj). 

Interpreting the effects of such a potential on the vortices themselves will lead us to 
studying a somewhat different problem. The dynamics of vortices can be extremely coun- 
terintuitive. Even in the simplest case of vortex lines — which are the only possible vortices 
in superfluids, but exist in ordinary fluids as well — the equations of motion for the lines' 
positions are first order in time-derivatives. This means that concepts like that of 'force' on 
a vortex (which we would naively derive from our potential) do not really apply. We remedy 
this by constructing an action that describes the vortex lines dynamics in the purely incom- 
pressible limit (sect. |6]). This action includes the aforementioned "dragging" interactions 
as long-distance (that is, non-local) Lagrangian terms. Then, our sound-mediated potential 
energy should just be interpreted as the first-order — in departures from the incompressible 
limit — correction to this action, whose consequences at the equation of motion level can be 
derived just by varying the action in the usual way. 

Finally, bringing everything together (sect. [?]), we rewrite the non-local dragging inter- 
actions as being mediated by an auxiliary local field, and we reintroduce the local couplings 
of vortices to sound, thus ending up with a convenient, local field theory describing the 
dynamics of vortex lines and their interactions — among themselves and with sound — to all 
orders in perturbation theory. We feel this to be substantial technical as well as conceptual 
improvement over the more standard "vortex filament" model [9j. 
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2 Coupling of incompressible vortex flows to sound 
waves 

We begin with the action for a relativistic perfect fluid. This can be written in two equivalent, 
alternative forms, depending on whether one treats, at any given time t, the fluid elements' 
physical ('Eulerian') coordinates x as functions of their comoving ('Lagrangian') coordinates 

0' [n], 

f=f(0^t), 1 = 1,2,3, (3) 

or vice versa [I2l [131 [1] > 

0' = 0'(x,t). (4) 

Even though the former viewpoint is more intuitive from a mechanical viewpoint, the latter 
makes the construction of the action more transparent, since the 0^s transform simply as 
three scalar fields under the spacetime symmetries (Poincare) acting on their argument 
x'^ = {t,x). There is also a number of internal (i.e., coordinate independent) symmetries 
acting on the (p^ fields, most notably the volume-preserving diffeomorphisms fi3[ [1] 



(5) 



corresponding physically to the absence of transverse stresses in a fluid. 

Following standard effective field theory logic [13], one finds that at long wavelengths 
and low frequencies, the action for the (j)^ fields is given by 



S = jd^xF{b), with b= ^detdf,(l)Wf'(l)-^ , (6) 

and where F is a generic function. A particular functional form of F corresponds to a 
particular equation of state. For instance, various thermodynamic quantities such at the 
energy density, the pressure, the temperature, and the entropy density, can be expressed in 
terms of F and its derivatives 



p = -F, p = F-Fbb, T = -Fb, s = b, (7) 

where the subscript b denotes differentiation. This formalism can be straightforwardly ex- 
tended to fluids at finite chemical potential for a conserved charge [H]. Notice that for an 
homogeneous, static fluid at equilibrium (at some given reference pressure), the physical 
coordinates coincide with the comoving ones, 

(0')eq. = ■ (8) 

This configuration has b = 1, which will be our reference value for b for what follows. 

Our goal here is to systematically incorporate compressional effects by expanding this 
Lagrangian around slow (with respect to the speed of sound) vorticose background fluid 
flows. At lowest order, this expansion was performed in [1], but we reproduce it here for 
clarity. The vortex/compressional mode separation is made most clear by working first in 
the x{(f), t) parameterization of the fluid, i.e. working in comoving coordinates. Additionally, 
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we make the speed of hght c exphcit in order to better keep track of relativistic effects. In 
these coordinates the action reads: 

S = -WqC^ J d^(f)dt det J/((det J-^)^/l-vyc^) . (9) 
Where J is the Jacobian matrix 

V is the fluid's velocity field, 

W = , (11) 

and wqc^ is its enthalpy density at the reference point 6 = 1: 

wqc^ = {p + p)b=i = -Fk{l) . (12) 

We are giving wq units of a mass density, and / is a dimensionless function of its dimensionless 
argument, defined by 

F{h) = -w,^f{h). (13) 

The above expression for the action is derived in Appendix |A] and was originally given in 
reference [13]. The form of equation (|9| is convenient for taking the non-relativistic limit. 
However, we will find that keeping track of all relativistic effects is not particularly difficult 
and so we will do so. 

For incompressible flows we have 

x(0,t) =xo(0,t) , (14) 

where, at each moment in time, XQ{(j),t) is a volume preserving diffeomorphism of the co- 
moving coordinates 0, 

det(Jo) = l, ^o^- = detf|§^ . (15) 



This flow will only be a solution to the equations of motion, given by varying (|9| with 
respect to x^{(j),t), when the time-dependence of the flow vanishes. Consequently, for slow 
(relative to the speed of sound) vorticose flows, xl^{(f), t) is only an approximate solution with 
compressional corrections. So for these nearly incompressible flows we parametrize the fluid's 
configuration by 

f(0,t) =xo(0,t)+^(0,t) (16) 

— * 

where ■ip{(f),t) vanishes in the limit of negligible time-dependence of a;o(0, t), and is curl-free 
with respect to Xq, in the sense that 

^f=V,,^{xo,t), (17) 

for some function \E'. It is easy to convince oneself that, at least to lowest order, this is 
the correct characterization of compressional modes. For instance the determinant of the 
Jacobian (10), which characterizes the compression level of the fluid, at linear order in ip is 

det J~ 1 + Va;o ■ V^, (18) 



5 



where we made use of the volume preserving property of Xo(0, t). So, a curl component (wrt 
to Xo) in ip does not contribute to compression — only a gradient component does. In the 
following we will denote derivatives wrt to xq simply by do — not to be confused with the 
usual relativistic notation of denoting do as the partial derivative with respect to time, which 
will be denoted here simply by df. 

We can now expand (|9| in terms of these small perturbations. To quadratic order in the 
compressional modes we have 

det J = 1 + [do^] + I {[do^]^ - [do^do^]) + . . . , (19) 

where [. . . ] represents the trace Meanwhile, the velocity field is given simply by 

v = vo + dttf{(P,t), (20) 

where vo is the divergence-free velocity field associated with xo{(j),t): 

vo = dtXo{(j), t) , do-vo = . (21) 

Taylor-expanding and changing coordinates from the (p to Xq, and dropping total-derivative 
terms, the action can be written to lowest order as 

S = So + woJ d^xodt {Ki^)' - Iclidoi^]^ + ^0 ■ iV^- livlido^P] + ...] (22) 



So = Wo / d^^dt ItTo 2 + . . . (23) 



where 



has no ^/'-dependence. The . . . 's in 5*0 are relativistic corrections, which are only powers 
of vo"^, c and Cg. Since So does not depend on ip, it should be thought of as the action 
for the incompressional flow parameterized by xo{(p,t). To use 5*0 in this sense, one should 
supplement it with a constraint that xq be a volume-preserving diff of (p. For instance, it is 
straightforward to check that adding to 5*0 the term 

j X{(P,t)[det{dxyd(f>i)-l] , (24) 

where A is a Lagrange multiplier, yields the correct Euler equation for incompressible fluids |^ 
However, as we stressed in the Introduction, the incompressional dynamics are very difficult 
to solve for, whereas, at least for slow fluid flows, we can make substantial progress by sys- 
tematically applying perturbation theory to the compressional degrees of freedom, for a given 
background incompressible flow. So, in the following we will treat the volume-preserving evo- 
lution xo{4>,t) as given, and we will solve for the small compressional perturbations of such 
a background. In the quantum field theory /functional integral language, this is equivalent 
to performing the path integral — at tree level — over the field, in the presence of given 
external sources xo{4>,t). 



^So [dip] = dilj)'- and [dipdil)] = djip'^diip^ etc. 

■^In the Euler equation, A plays the role of pressure. 



6 



In (22), the -^'s remind us that the partial time derivatives were taken at fixed comoving 
coordinates, i.e. fixed (p. These are the usual convective derivatives of fiuid dynamics. In the 
Xo coordinates, they simply act as 



D_ 
Dt 



g{xo, t) = dtg{xo, t) + {vq ■ do)g{xo, t) . (25) 



This expression is exact, to all orders m ip: it is a mathematical identity, which follows 
simply from the chain rule; the fact that Xo{(p,t) only reproduces part of the full, physical 
fiow, is irrelevant. Finally, we have used that 

/'(I) = 1 , r(l) = c'Jc' (26) 



(cs is the sound speed), which follow, respectively, from the definitions (12), ( |13|), and from 
a straightforward analysis of the perturbation spectrum about the background (|8j) IT]. 

And so, keeping the free quadratic terms in ip (those not coupled to vq), and the linear 
coupling of ip to the background vortex velocity, we finally write the lowest order action 
containing the compressional modes as 

AS = Sfrcc + 5'int (27) 
Stee = WoJ d'xodt \ ((9*^0' ' c'[5o^]') (28) 

5int = wa j d^xdt{vl{vo ■ do)tpi - l§vl[do^] + ...). (29) 

We have discarded a vq ■ dtip term after noting that, as ip = do"^, we can place this spacial 
derivative on vq by a simple integration by parts, and the term will vanish as 9o ■ = 0. 
This simple action is enough to generate the famous "analogy" of Lighthill [1], with an 



additional relativistic correction. When we vary (27) with respect to ip'' we recover Lighthill's 
limit of (|2|). Of course, for the fiuids Lighthill was implicitly considering the relativistic effects 
would be extremely subleading, however for ultra-relativistic fiuids where can be as large 
as this correction term can become a leading-order effect. 



3 Sound emitted by a vorticose source 

As an illustration of the utility of the effective field theory techniques, we reproduce Lighthill's 
result with additional relativistic correction terms. Instead of manipulating the lowest order 



equations of motion we can work directly from the action given by (27) utilizing the stan- 
dard tools of quantum field theory to calculate physical observables. In particular, using the 
language of Feynman diagrams and restricting to tree level amplitudes allows us to isolate 
classical quantities [21] . For comparison, a derivation of Lighthill's result using more stan- 
dard techniques can also be found in Landau and Lifshitz's classic Fluid Mechanics textbook 

m 

Say that we have a vorticose, possibly turbulent, configuration v{x, t) with typical fiuid 
fiow oi V <^ Cs and some characteristic size I. The velocity field acts as a source for sound 

^ It should be noted that there is a small typographical error in their final result. 
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waves. We can formally deal with this emission process via QFT tools, by computing the 
amplitude for emission of a single "phonon" from our compact vorticose source. The linear 




(in ^p) terms of Siat in (27) contribute to the "tadpole" diagram 

= tM{p), (30) 

where A4.{'p) denotes the probability amplitude to emit a single phonon of momentum p. 
Given the probability amplitude one can calculate the emission rate with the standard for- 
mula: 

where T — i- oo represents the total integration time, and will drop out in time averaged 
quantities (which is what we are after). The standard relativistic formula has also been 
modified with c — )■ as the energy of a phonon with momentum p is simply CsP (see ref. [T] 
and Appendix [b| for details). 

Proceeding in the usual manner one finds: 

^M{p) = w'J'^k*^{csP,p) (32) 

where Hij {x, t) is the "kinetic tensor" 

Kij{x, t) = ViVj - \5ij ^ , (33) 

and its Fourier transform, 

kij {uj,p)= [ £xdt (f , t) e'^^'e-'P-"^ . (34) 



Before inputting this amplitude into equation (31), we can make use of the multipole ex- 
pansion: The typical frequency of the vorticose source is w ~ The emitted sound waves 
inherit such a frequency, and have therefore a typical wavelength A ~ Cg/u ~ i ■ [cg/v), 
which is, in our approximation, much bigger than the vorticose source's size i. This means 
that we can treat the source as point-like, or the sound waves' wavelength as infinite. In the 
monopole approximation we have 

kij{csp,p) ~ kij{csp,0) . (35) 
In this approximation, integrating over the solid angle we have 

U^du Wo r .2 , o\~ |2l 

aL[uj) = — — r Kjj + 2\Kii\ , (3d) 

where the k's are evaluated at frequency u and vanishing p, and we have used n^n^ = ^5ik 
and riinknirim = j^iSikSim + SuSkm + ^im^ki) (the 'overline' denotes the average over the solid 
angle.) 
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The emission rate is not the most natural quantity to consider when calculating classical 
wave emission, a more typical one would be the power radiated. It is easy enough to augment 



(36) in order to get this more standard quantity. If we have F = J dT we can write the power 
as simply P = j E-dT. In our case, the energy of a single phonon is simply u (we are working 
in h = 1 units). Extending for convenience the range of integration to negative frequencies 
(and dividing by 2), the total radiated power is 

(37) 



T 1207r2c! 



I ^ii I ~l~ ^ I ^ij I 



Noting that, for real /(t) and g{t), 

we see that the inclusion of the 1/T factor, where T — )• oo, gives us time averages. Meanwhile, 
the restriction to vanishing p gives us volume integrals of the form 

Kij{t) = j d^xKij{x,t) . (39) 

Putting everything together, for the radiated power we finally get 

Wo 



P 



(^^) + 2(K./) . (40) 



where (...) denotes the time average, and the over-dots time derivatives. 

To compare this result to the standard one — see e.g. [2] — , we define the traceless and 
pure-trace quantities 



Qij{t) = j d^x {viVj - Iv^Sij) , Q{t) = j d^x , (41) 
which we can use a basis for our tensor Kij. In terms of these we get 

'l{Q'){l-4 + -J) + UQ^/)] ■ (42) 



p= ^ 



This matches the result contained [2] with the addition of relativistic corrections. Note that 
these are much smaller than the leading terms only if Cg is much smaller than c, regardless 
of how small v/c is. For relativistic equations of state with ~ 1/3 c^, the relativistic 
corrections are unsuppressed. The above expression is correct to lowest order in v/cs — 
which is our small expansion parameter, associated with near incompressibility — but to all 
orders in c^/c, and thus applies to those cases as well. 



By going to higher orders in the diagrammatic expansion of eq. (30), one could system- 
atically compute higher order (in v/cg) corrections to this result — for instance "radiative" 
corrections, like was done in the case of gravitational wave emission by binary systems in 
[16] . The advantage of the EFT techniques for performing these higher-order computations 
is manifest: for instance, the UV divergences associated with the point-like approximation 
for the source can be handled in the standard and well understood ways of renormalization 
theory. 
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Figure 1: Feynman diagrams contributing to the scattering of sound by a background fluid flow. 
The circled crosses denote tp-vo interaction terms. There is no propagator associated with the double 
lines, which just depict the external source vq. 



4 Scattering sound waves off vorticose sources 



If we include higher order terms in (29) we can calculate other interesting quantities. For 
example, the scattering of sound waves by a nearly incompressible flow. As emphasized in 
[18] this can be a powerful probe of the details of the fluid flow, e.g. to address questions 
about turbulence. The simplest scattering diagram will come from terms in the action 
quadratic in — one ip for the incoming wave, one for the outgoing one — and coupled to vq. 
In principle there are other diagrams contributing to the same physical process, depicted 
in fig. [T] However, it is easy to convince oneself by simple power counting that, for small 
vq/cs, the quadratic-in--?/^ vertex is the leading contribution: as derived in Appendix |A| it 
involves a single power of vq, whereas all the other diagrams start at order Vq. The 'cross' 
in the second diagram of the figure yields directly a Vq factor, while the crosses in the other 
diagrams yield a factor of Vq each, but there are two such vertices per diagram. 



Continuing the expansion that led to (29 ), and which is summarized in detail in Appendix 



[Af we obtain for this interaction term in the action: 

Sint DS^,n=wo [ Sxdt \{dtiJi){v ■ ViP') - '^-^(vA'ip'Wip] + 0{v') \ . (43) 



C2 



The second term is a relativistic correction, which is important if the speed of sound is not 
much smaller than that of light. With an abuse of notation, we have stopped differentiating 
between x and xq, and between v and vq: all x's and w's should be thought of as being xq's 
and fo's. 

When calculating the power emitted by a vorticose source in the previous section we only 
made one assumption about the flow (in addition to its having a finite size): its flow velocity 
was very subsonic, with v/cg <^ 1. Now however, we have another dial we can control, that of 
the frequency of the injected sound waves. In particular, we take the same limit as [18] and 
take the frequency of the incoming sound to be much larger than the natural frequency of the 
vorticose fluid flow: uj^^ u^. This is an experimentally convenient limit, as it distinguishes 
the scattered radiation from the naturally emitted one discussed in the previous section. 

For given velocity profile ^(x, t) of the source, the (co'i,^!) — )■ (a;2,P2) scattering amplitude 
associated with the interaction terms above is simply 

iM = -i| (pi ■ P2) [(^ipi + UJ2P\] - ^ [(^ip\ P2 + Pi] '^v\Auj, Ap) , (44) 
where Au and Ap are the energy and momentum transfers: 

Aid = cji - CU2 , Ap = pi - p2 . (45) 
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According to our assumptions, the typical frequency and wave-number of f * are 



Au ~ -C 



Ap 



(46) 



where i is the typical size of the source. We thus see that from the sound waves' viewpoint, 
the energy transfer is negligible, both in absolute terms and relative to the momentum 
transfer: 

Auj Au V Ap Ap 

<1, ~ — . 47 

Ui Ui Cs Pi Pi 



To simplify the amplitude we can thus set 



|Pl,2| ^Uj/Cs . 



(4J 



On the other hand, the momentum transfer — and thus the scattering angle — can be sizable, 
as long as oj^/uj is not much smaller than f/c^ : 



Pi 



V U 



(49) 



We are thus left with 



iM 



iPi • P2) - % (Pi + P2) ■ v{Auj, Ap) . 



(50) 



Inserting this into the standard formulae for the scattering cross section, which, for the 
convenience of the reader, are also quickly derived in Appendix [B| we have 



da 



dVLd{Auj) 4c6 (27r)3 T L 



{pi ■ P2Y - 7# (Pi ■ P2) + 



\{Pi +P2) ■ 



(51) 



where T — t- 00 is a long observation time — long enough to observe the whole scattering 
process (the same T should be used to define the Fourier transform v). 

If the source's velocity field is nearly stationary, v ~ (2-7?) 6{Auj), or more generally if one 
does not have the experimental resolution necessary to detect the small frequency transfer 
Au, one can integrate in Ao; and end up with a time-averaged expression: 



da 

dn 



1 

4rf 



(2vr) 



iPl-P2] 



2ci 

r2 



iPl ■P2) + $ iPl + P2)\PI + P2y{v"v^) 



(52) 



where now v = v{Ap,t) denotes a purely spatial Fourier transform of the source's velocity 
field, and (...) denotes a time-average. Needless to say, the second and third terms in brackets 
are relativistic corrections, and — like in the previous section — they are exact to all orders in 
Cs/c. 

The result expressed in (51) can be recast in a form that can be compared to previous 
results obtained by more traditional (and much more laborious) methods [18J. However, 



the comparison takes some non-trivial translating and so it is contained in Appendix C.l 
When the smoke clears, we see that our answer disagrees with that given in 



by an overall 

factor of 2. To test the self-consistency of our computations, we check whether they obey 
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the optical theorem. Since the optical theorem relates something quadratic in the scattering 
amplitude (the total cross-section) to something linear in it (its imaginary part in the forward 
limit), this is a non-trivial check of the overall normalization of our final result. According to 
standard cutting rules, at this order in perturbation theory the optical theorem should relate 
the cross section we have just computed to the imaginary part of the the third diagram in 
fig. [T| As we show in Appendix C.2[ the check works out, which gives us confidence in our 
results. 



5 Sound- mediated long distance interactions between 
vortices 

In addition to the physical processes considered above, which have been previously treated 
by perturbation theory performed at the level of the equations of motion (at least to lowest 
order for non-relativistic fluids), our formalism naturally motivates the calculation of the 
long-range interaction due to the exchange of compressional modes between two physically 
separated vortex configurations. The effective field theory approach invites such a question, 
and gives it a clear-cut qualitative answer: our vortices interact with compressional modes; 
these are gapless, and, as a consequence, can be exchanged between arbitrarily distant vor- 
tices; at some order in perturbation theory, this exchange will yield a non-trivial long range 
interaction between vortices. There is no question about it. For instance, this is how the 
classical 1/r Coulomb interaction arises in QED — via the tree-level exchange of virtual pho- 
tons. On the other hand, such a question is not a natural one to ask in the standard approach 
to fluid dynamics, and in fact, to our knowledge, these long-range interactions have never 
been postulated before — let alone computed. 

Notice that there are other long-range interactions between vortices, which survive even 
in the strict incompressible limit [cg — ?■ oo), and are well understood in purely kinematical 
terms. As we will describe at some length below, each vortex typically carries a long distance 
"tail" in its velocity profile, scaling for instance as 1/r for infinitely long vortices. Then, if 
one is given several well-separated vortices, each will be dragged by all the others' velocity 
tails. As will see, our sound-mediated interaction will contribute a small correction to this, 
suppressed by (f/c^)^, but by no extra powers of r. 

In [1] we have initiated the analysis of our long-range interactions, by solving the lin- 
earized equation of motion for the sound modes and plugging the solution back into the 
action. Here we will use a more systematic technique — that of the effective action]^ Even 
though at the classical level using the effective action is equivalent to solving the (generally 
non-linear) equations of motion for some of the fields and plugging their solutions back into 
the action, and even though in the end we will restrict to a lowest order computation thus 
reproducing the result of [1] , we find it useful to set up a more systematic framework anyway, 
to pave the way for higher-order computations (which will appear elsewhere [23]) or for more 
general ones. 

For our system, the sources will be given by vortices characterized for simplicity by the 
^ An excellent pedagogical introduction to the effective action as utilized in a very similar context is given 

by m- 
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same typical length scale / and the same typical velocity ^ c^, separated by a much 
bigger length scale r ^ Z, so that they appear as point-like sources to each other. These 
vortices interact via the exchange of compressional modes — virtual phonons, in our QFT- 
inspired language. We can compute the potential energy due to this exchange, and standard 
thinking indicates that its gradient will give the "force" that the sources exert on each other. 
As we will see, this intuition is fundamentally mcorrect, but the potential energy is still 
a well-defined physical quantity worth computing, and we can worry about its dynamical 
implications later. 

The relevant quantity is the effective action Ses one is left with after integrating out the 
compressional modes in the path integral. Schematically 

Since is gapless, S^s will contain non-local interaction terms — precisely what we are after. 



As it turns out, the simple lowest order action given by (27), 

S[Xo,1p] ~ So[Xo] + Sireel'ip] + Sint[Xo,ij] , (54) 

is all we need to perform this calculation to leading order. 5*0 is the lowest order action 



for xo(0, t), given in (23). S^ee is the quadratic action for ip, and ^int contains its linear 
couplings to Vq: 



Sint = J d^xiP-J , J{x) =wo{-vo-{vo- V)vo + l"^V{vo^)) ■ (55) 
Therefore, to lowest order we just have 

Seffixo] = So -l\nZ[J] , Z[J] = J l)^^e*5fr,,+i/d4,fj+e terms ^ (gg) 

and we can compute the Gaussian functional integral via the standard completion of the 
square method. Up to an irrelevant constant we get simply 

Sos = So + l [ d'^xd'y Ji{x)A'^{x - y)J,{y) , (57) 



2 , 

where A*-' is the usual Feynman propagator defined by 

AH.^y).^{0\T{Ay)r'(.)} |0> ^ - / p^z^j^^ . (58) 

Diagrammatically, this leading order correction to the effective action is given by the simple 
diagram denoted in fig. [2j 

The compressional modes that mediate a long-distance interaction between slowly-moving 
sources are going to be (very) off-shell, i.e. — /cq +Cgfc^ ^ 0, and so we can take e — )■ straight 
away. In particular, we expect the relevant frequencies and momenta to be |^ 

~ y, ~ - , (59) 



^We are assuming that v v'^/l, so that the typical frequency of the sources is v/l. However, in many 
cases, Hke for instance those we wiU discuss below, one has nearly stationary, incompressible flows, in which 
case the typical frequencies are zero in first approximation, and the associated "virtual" compressional modes 
are maximally off-shell. 
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Figure 2: Lowest order diagram contributing to the potential energy mediated by virtual sound 
between two vorticose sources. 



and so if we take v/cs ^ l/r we can expand the denominator of the propagator as a power 
series in [k^/cgk): 

' ' ^1 + ^ + ...). (60) 



-k^ + Clk^ Clk^ V 

This is analogous to the low-energy expansion of a massive propagator in relativistic QFTs, 
and for us it means that our interactions, despite being genuinely non-local in space, can be 
expressed as a series of interaction terms that are at least local in time. 

So, to lowest order in this expansion and isolating the interesting interaction term (that 
is, the cross term when we write the total vortex velocity as the sum of the two separate 
sources v = Vi + V2) we have simply, 

S,^ = So + AS = So + - d'xdS^J\{x)Ji{y)^^e^^^--y^ (61) 

Wo J (27r)4 cffc4 

= + ^° (0^ ^^^^'^^"^^ ' (62) 

where we have freely integrated by parts and utilized the divergence free nature of and 



the k's are the two sources' kinetic tensors, as defined in (33). 

Now, as these sources are supposed to be localized (that is, their velocity field falls off 
sufficiently fast away from their respective centers) we can multipole expand each source, 
keep the lowest order term in {l/r) (the monopole), and perform the k^ integral, which 
generates a (27r)5(ti — ^2)- We may therefore write the integral of interest as: 

A5 ^ ^ / rft Kl\t)Ki\t) [ 4^ 'l^^e^'^-r ^ (63) 



where r is the vortices' relative position vector, and, like in sect. |3j each K is defined simply 
as the monopole moment of the corresponding k, 

K'^ = j d^x k'^ = j d^x {v'v^ - \%b'^ v^) . (64) 

The k integral is straightforward to perform by first rewriting the fc's at the numerator as 
gradients w.r.t. r, and then noticing that 



/ 



(^k \ ^ 1 

g«fc r _ _ ^ _|_ iR-divergent, f-independent pieces, (65) 



(27r)3 k^ Stt 
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as can be checked by introducing any IR-regulator, like e.g. — i- 1/(A;2 + m^)^. 

Recalling that the potential enters the action with an overall minus sign, we can finally 
write the effective vortex-vortex potential mediated by sound in a compact form as 



A V = K'^Ki' didjdkdi r . (66) 
The derivatives, once expanded, yield 

didjdkdi = \ [— {^ij^ki + 2 perms.) — 15 ■ f^f^f^f'' + ?>{r'^f^ki + 5 perms.)] . (67) 



Notice that the monopole moment (64) obeys positivity properties that forbids its van- 
ishing as soon as one has some v — no matter how inventive one is in devising such a velocity 
profile. For instance, its trace is 



K'' = I d'xv'{l-l§) , (68) 



which, using c^/c^ < 1/3, is always bigger than f ^v"^ — the kinetic energy per unit enthalpy. 
As a result, the monopole- monopole interaction that we have computed will always be the 
most important long-distance interaction mediated by sound 

For more general or higher order computations it is probably more convenient to perform 
the multipole expansion directly at the level of the action, before integrating out tp, along the 
lines of [8]. In this way, on top of organizing the perturbative expansion in a more systematic 
fashion, one is able to handle the UV divergences that unavoidably will show up at some order 
in perturbation theory using the standard tools of renormalization theory. For instance, 'self- 
energy' diagrams — diagrams in which a vortex exchanges a sound mode with itself — are UV 
divergent. Their divergent contribution to Ses can be reabsorbed via renormalization into 
the coefficients of terms that are already present in the (multipole-expanded) Lagrangian. 
We leave a systematic treatment of the multipole expansion for future work. 



Eq. (66) is our final result for the long-distance, monopole-monopole interaction be- 
tween localized vortex configurations. It is easy to estimate its size. The individual kinetic 
monopole moments are roughly K"^^ ~ v'^l^, the derivatives structure acting on r scales as 
r~^, and so the effective potential scales like 

Ay ~ ^ ^ ■ E^n , (69) 

where i^kin ~ it the typical kinetic energy of the vortex configuration. We will see 

in sect. |6]that the aforementioned kinematical dragging phenomenon can be modeled via a 
long-range potential energy that, for localized vortices, scales as (//r)^ ■ -Ekin- Our effect is 
thus suppressed with respect to this one by an extra {v/csY factor. 

5.1 Potential between two vortex rings 



In order to get a better physical sense of eq. (66), we can compute it for the simple configu- 



ration of two interacting circular vortex loops — the so-called vortex rings. Vortex rings are 
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not only the simplest localized vortex-line configuration; they are also beautiful, fascinating 
objects of intense experimental interest and study, for humans P,|2l] as well as other species 
[25] (see fig. [s]). For the following we will ignore the relativistic correction term in K"^^ , but 
its inclusion is straightforward. We refer the reader to [9] for an extensive review of the 
vortex rings' physical properties. 

Say we have a vortex ring lying in some plane with normal vector h. The radius of the 
ring is given by the circulation by F, and the vortex line thickness by some cutoff a (the 
so-called core radius). By symmetry K"^^ will take the form: 

K'^ = A5'^ + Bh'n^ . (70) 

The parameters A and B are given by the scalar integrals 

^ = \ j d'x{v^ -{v-hf) (71) 
B = lJ d'x{iv-hf-y) (72) 

For distances much larger than R the velocity flow is given by that of a dipole, that is, 
the velocity falls off like 1/r^, and so, as anticipated and desired, the main contribution to 
K'^^ comes from the region near the vortex ring itself. However, at very close distances to 
the vortex line the same integral will diverge logarithmically (if the core radius a is taken 
to zero). For a ^ -R, we can thus get a good estimate of the integrals above by isolating 
the coefficient of log a, which we can compute in the straight vortex line limit. Then, by 
dimensional analysis, log a has to be accompanied by — log R. Geometric factors coming 
from precisely integrating over scales of order R (which we can do, at least numerically) and 
from the details of how the UV cutoff is actually implemented (which is model-dependent), 
will show up as finite, order-one additions to this large universal log R/a contribution, and 
are therefore sub-leading. Implementing this trick we find: 

K ~ IT^R \og{R/a) {I + n ® n) . (73) 

However for realistic vortex rings in normal fluids (as opposed to superfluids), R is not that 
much bigger than a, and the "divergent" log is so weak that one cannot neglect the allegedly 
sub-leading finite pieces. Therefore, in order to be completely general we can write 

K = T^R{aI + Pn(^n) . (74) 

where a and P are order-one (or log) coefficients that in general depend on the exact core 
structure of the vortex rings. 



Inserting this into (66) and assuming for simplicity the same core structure for the two 



rings, we can write the effective potential as 

Al^ = £^ iTlR^) (TlR^) ■ ^ ■ /(n^, n^, f) , (75) 

where / parameterizes the angular structure: 

/(ni, n2, f) = -(/?' + 4a/3) - 2/3\n, ■ h^Y - 15P\h, ■ ffih^ ■ (76) 
+ (6a/3 + 3/3^) [(ni ■ ff + (na • ff] + 12(3\ni ■ r){n2 ■ r){ni ■ ns) . 
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Figure 3: Dolphins can set up 
very clean vortex ring configura- 
tions, and play with them. They 
blow air into them, which, given 
the pressure gradients associated 
with the fluid flow, 'floats" to 
the core of the ring, thus serving 
beautifully as a tracer. Pictures 
taken from W, 



Notice that / can be either positive or negative. For instance, if we take a = /3 = 1 for 
simphcity, then for a configuration with fii = 7x2 = r one gets / = 8, while for hi = n2 ^ r 
one gets / = —7. That is, depending on the geometric configuration, this potential can be 
both "attractive" and "repulsive". As we will see in a moment, the quotes are in order, 
because these terms are extremely misleading for vortex dynamics. 

5.2 How (not) to interpret this potential 



For standard systems, one usually thinks of a potential like that given by (66 ) as a function of 



r'that will generate a force in the usual way, F = The problem with this interpretation 

for us is that the notion of "force" does not really apply to vortex dynamics. To see why this 
is the case, let us restrict ourselves to the study of vortex lines, which exist both in normal 
fluids and in superfiuids, and of which the vortex rings we just discussed are an example. 

The way that we visualize a vortex line is by thinking of it as exactly that, a line, or 
more generally a curve. But really such a curve is a placeholder for a very special extended 
field configuration: the curve is the locus where the vorticity is nonzero, but the velocity 
field extends to large distances from the curve, typically in a 1/r fashion for very long lines. 
As a consequence, you cannot do to a vortex line the usual kinds of things that you would 
to a string-like object in empty space, like, for instance, boost it. Doing so will change the 
boundary conditions at infinity of the velocity field, thus in effect boosting the whole fluid. 
This does not mean that vortex lines cannot move, however, it just simply means that for 
the velocity field to go to zero at infinity, the vortex line will move in only one particular 
way; there is no free initial condition for their velocity so that, at any given time, the local 
velocity of any vortex-line infinitesimal element is completely determined by the dynamics. 
For instance, isolated, perfect vortex rings always move at a constant speed determined by 
their radius and circulation (see for instance [H]). Why do they do this? 

In the incompressible limit, as p is no longer a dynamical variable, the hydrodynamical 
equations of motion are for v only and are given by the simplified continuity equation and 
the curl of the Euler equation: 

V-v = (77) 
d 

— cJ = V X (vxu) (78) 
dt ^ ^ ^ ^ 



where cJ = V x v is the vorticity field. The linearity of (77) implies that if we have multiple 
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isolated vortex lines (lines where the vorticity is non-zero), then the velocity field outside 
these lines obeys linear superposition. Furthermore, Kelvin's theorem (which follows from 



(78)) applied to infinitesimal comoving loops surrounding the vortex lines, implies that each 
line will move with the fiuid fiow generated by itself and every other vortex line [?] . Since the 
relation between cJ and v is formally identical to that between current density and magnetic 
field in magnetostatics, and since for vortex lines u is given by a sum of circulations times 
delta functions supported on the lines, with the circulations playing the role of total electric 
currents fiowing along these lines, the local velocity of the fiow at x generated by all the 
vortex lines can be written in a Biot-Savart fashion as 

m = -y.^ [ fcS X dx' (79) 



^ 47r J„ \x - x'P 

where r„ is the circulation associated with the n-th vortex line 

As we can see, given that a vortex line element placed at x will move with this velocity, 
the vortex lines' positions obey first order equations of motion, in contrast to the usual 
second order equations of motions of any string-like object in empty space. This is the root 
of their peculiarity: with first-order equations of motion, there is no room for "forces" , and 
the effective potential we have computed above must then be interpreted more carefully. 
Our strategy in the next section will be to derive a Lagrangian formulation that reproduces 
precisely these first-order equations motion. Then, our sound-mediated effective potential 
will just be a correction to such a Lagrangian, which will yield a corresponding correction 
to the equations of motion via the variational principle. 

The model that we briefly sketched above is called the vortex-fllament model and was 
pioneered in [22] and [30] • A modern review of such topics is discussed nicely in [S]. In the 
next sections we will flnd a reformulation of the same model that we believe is likely to lead 
to substantial progress in solving systems involving vortex lines. 



Before we proceed, note that the integral over the "self- interacting" line in (79 ), is going to 
diverge logarithmically as x' — )■ x. Introducing a flnite UV cutoff a, say the line's thickness, 
one gets that the leading contribution to the self-interacting integral is 

^ ^ X dx' z^\og{R/a) + ... , (80) 



self 



Ix — x'P R 



where R is the line's local radius of curvature at x, z is the direction orthogonal to the plane 
deflned by the local curvature, and the dots stand for subleading pieces, which come from the 
integral over regions that are far from x, and which are flnite for a — )■ 0. This UV divergence 
tells us that there is a mild dependence of the velocity of vortex lines on the details of their 
core structure. Changing these details, or equivalently changing the deflnition of a, will only 
affect the subleading, flnite pieces. As usual, the coefficient of the log divergence is universal. 
If one is in the large \og{R/a) limit, the velocity is dominated by this universal piece. In 
the vortex fllament model this approximation is called the local induction approximation, 
because the local velocity is mostly induced by the local curvature of the vortex line. 



''There is a well-defined circulation r„ = jg oj-dS for each line (the integral is taken over the cross-section 
of the "line"), which is in particular constant along that line, because of Stokes theorem, and constant in 
time, because of Kelvin's. 



18 



6 An action for vortex-lines and their interactions 



Given that this paper focuses on the Lagrangian formahsm, it is natural to ask if there 



is an action that will generate (79) upon variation. As mentioned above, a byproduct of 



finding such an action is that it would provide us with a straightforward interpretation of 



our sound- mediated potential (66): this should simply be interpreted as a correction to the 



action, with obvious implications for the equations of motion. 



Strictly speaking, the lowest order incompressible action (23), Sq = wo J ^vq"^, contains 



the dynamics we want to reproduce. However, we want to restrict the velocity field to 
configurations with a vorticity field localized on lines, and we would like to use the positions 
of the lines themselves as degrees of freedom. It is not trivial to rewrite 5*0 as an action for 
these degrees of freedom. So, for the moment we will just guess what the right action is, and 
we will comment later on what its relation to 5*0 is. 

Consider for the moment just one vortex line, parameterized by X{t, A), where A is some 
parameter running along the line, over which we are going to integrate. The action that we 



are going to write down should produce the equation of motion (79) upon variation w.r.t. to 
X. Let's focus on the left hand side first. To get the velocity, which is just v = dtX, one 
needs a term in the action with one time derivative only. A nontrivial, natural candidate in 
2D (for vortex dots) would be 



J dt liX^dtX^ - X^dtX^) = j 



dt\e'^ X'dtX^ 



Varying w.r.t. to X^ yields dtX"^, and varying w.r.t. to yields —dtX^. The antisym- 
metric structure is crucial in order not to end up with a total derivative. Then, a natural 
generalization to 3D seems to be 



dtdXy^^X'dtX^dxX^ 



(82) 



The vector d\X is tangent to the curve, and so the combination e^^^d\X^ behaves like the 
two-dimensional e-tensor for the plane locally orthogonal the curve, thus giving us back the 
2D expression above. The only apparent downside is that when we vary w.r.t. X to get the 
equations of motion, we now end up with an extra d\X^ multiplying (via a cross product) 
the velocity on the left hand side. 



Let's ignore this for the moment and let's move to the right hand side of (79). It contains 



a line integral over all vortex lines, including the one we are trying to understand the motion 
of. Then, the corresponding term in the action should be represented by a double line 
integral, over all possible pairs of lines, including those pairs made up of the same line taken 



twice. Notice that the line element in (79) can be written as dx' — >■ d\d\X. Notice also 



that eq. (79) involves a cross-product. Given the extra cross product with dxX we expect to 
get from (82), and given that, schematically ee ~ 5, we end up with the following inspired 



guess for the action of a system of vortex lines: 



S = C dt 



J2 I dXlT^^^'Xld.XidxX'^ + ^Yl I 

n n,m 



dXdX Lj^r^ 



dxXn ■ d\>X„ 
\ Y Y \ 



(83) 
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where the sums run over all the vortex lines, C is an overall — yet to be determined — constant 
(whose value does not affect the equations of motion), while the relative l/Svr is needed to 
get the correct equations of motion, as we will show in a moment. In this expression for the 
action it is understood that, in the double-integral term, X„ is parameterized by A and Xm 
by A'. As we stressed above, the double sum also includes = m, in which case the same 
line is traced out twice, independently by A and A'. 

Varying the curves in the usual way, X„(A,t) — )■ X„(A,t) + 5X„(A,t), and ignoring 
boundary terms, we get the equations of motion for the n-th line: 

e'^^dtXidxX'^ 

' (X„ - X^YdxXi dyXl (X„ - XmYdxXi dx^Xi 



r 

47r 



dX' 



IXr, — Xr, 



\X„ 



X |3 



4) 



or, in vector notation (using the identity dx{bxc) = b{d-c) — c{d ■ h) on the right hand 
side): 



dtXn X dxX„ 



d\ 



, {Xn — Xm) X dx'X„ 
\X - X |3 



X dxX„ 



i5) 



Now, it seems that our action rtSSft has not successfully reproduced the desired eom (79); 



there is an additional cross product with dxXn which we cannot simply remove from both 
sides of the equation. As a result, the equations of motion generated by our action are a 
little less stringent then the ones generated by "physical reasoning" . Should we be worried 
about this? 



As it turns out, the little extra freedom in (85) makes manifest a particular redundancy 



of how we describe the physical configurations of our system, which so far we have been 
reticent about. Given the spatial configuration of the vortex lines in our system at some 



particular time t, eq. (85) will dictate the velocity of a particular point on the vortex line 
only in the direction normal to that line. That is, the component of the velocity parallel to 
the curve is undetermined. But as the fluid velocity flow is completely determined by the 
positions of the vortex lines, any motion that does not change the locations and shapes of 
the lines, such as motion along the lines, is totally unphysical. In hindsight, it is obvious 



that this is related to a reparameterization invariance of the action (83): 

A ^ A(A,t) . 



(86) 



Both terms in the action are invariant, because the pieces that transform non-trivially are 
combined in structures that are manifestly invariant: at fixed t one has 

dX dxXn = dXn , dX dX' dxXn ■ dyXm = dXn ■ dXm . (87) 

This is a form of gauge invariance, which must be treated in the usual ways; for instance, 
in order to solve the equations of motion, one should first choose a gauge. The traditional 



form of the eom, (79), corresponds to a particular (and not so simple) gauge choice. Perhaps 
a more "physical" choice in our language would be to choose dtX ■ dxX = 0, i.e set the 
parallel component of the velocity to zero, for all time. But depending on the problem at 
hand, other choices might be considerably more convenient. For instance, to study an infinite 
vortex line that is roughly aligned with a coordinate axis, say the z axis, one can choose the 
gauge X = z = X^, thus effectively eliminating the X^ degree of freedom from the problem. 
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6.1 Fixing the overall coefficient 

While we have been able to reproduce the dynamics of the fluid system given the action 



(83), we still have to determine the action's overall scaling. We can fix this constant C by 
matching a physical quantity for a simple physical configuration. As usual for effective field 
theory, the matching of Lagrangian coefficients can be done via an observable of an idealized, 
simple configuration, and the inferred values for the coefficients can then be used for all other 
configurations as well. 

Consider a single vortex line aligned with the z-axis, with circulation F. Its velocity field 

is 

^ r 1 

v=—-(p. 88 
zvr r 



The original non-relativistic action for an incompressible fluid, eq. (23), implies that such a 
configuration has an energy per unit length 



2 



— = tf / lvl(f X = Wo— log L/a , (89) 

where L is some IR cutoff (e.g. the size of the container) and a some UV cutoff (e.g. the 
vortex-line thickness) |^ 



We can now compute the same physical quantity with our new action, eq. (83). The first 
piece — the 'kinetic' term — does not contribute to the Hamiltonian of the system, since it is 
linear in X: 

dL 

E = H = X' -L . (90) 

dX' ^ ' 

We are thus left with the second piece only, which yields 

E = -C— / dzdz- r , (91) 

87r y k-^'l 

corresponding to an energy per unit length 

dE 

- = -C-logL/a. (92) 

Notice that, given the different geometric nature of the integrals involved, the cutoffs used 
here cannot have exactly the same meaning as those used above; the universality of log 
divergences allows us to be cavalier about this. Comparing the two expressions for the 
energy per unit length we get simply 

C = -Wo . (93) 



^As before, as long as we are interested in the log-divergence only, we need not be specific about how 
these cutoffs are implemented, since the coefficient of the log is completely independent of these details — only 
additive finite pieces are sensitive to them. 



21 



7 The hydrophoton 



Eq. (83) reproduces the correct dynamics and energetics of generic vortex line configurations, 



and is thus the correct action for them. Nevertheless, it has an obvious annoying feature: 
due to the double A-integral in the second term, it is a non-local action. To develop an 
intuition about the (counter-intuitive) dynamics of the system, especially for perturbation 
theory questions (stability of solutions, properties and interactions of small perturbations, 
etc.), having a local action would be much more convenient. 

Fortunately, given the 1/r nature of the non-locality in question, we know exactly how 
to fix the problem. We 'integrate in' an auxiliary local field A{x, t), coupled to the 'currents' 



r„ dxXn, which are localized on the vortex lines: 



'S'local — Wo 



J2^n [ dtdXle'^'^Xi^dtXidxX^ 
+ I d^xdt i (diAj) ^ ~^^n I dtdX dxXn ■ A{Xn, t) 
One can solve the equations of motion for A deriving from this action. 



(94) 



V^A{x, t) 



dxrr.dxXnS'{x-Xn{x,t)) 



(95) 



and plug the solution back into action, thus getting back the non-local action (83). Since A 



appears at most quadratically in the action, this equivalence is exact, even at the quantum- 
mechanical level, which might be important for zero-temperature superfiuids. 

In this new local picture, interactions can be described in standard field theoretical terms: 



the action at a distance in (83 ) has been replaced by local interactions of each vortex-line with 



A. Then, different vortex-lines, or different pieces of the same line, interact by 'exchanging' 
A. It is obvious what A corresponds to in the magnetostatics analogy of the previous section: 
it is the vector potential in Coulomb gauge (V ■ A = 0). Indeed, given the (localized) current 
density 

J{x,t) = Y,^r^j dXdxX^6''{x~X^{X,t)) , (96) 



the j J ■ A term in (94) is precisely the standard interaction between current and vector po- 
tential, and eq. (95) is precisely the equation of motion for the vector potential in Coulomb 
gauge, V'^A = J. We dub A the 'hydrophoton' field. For more general vortex configura- 
tions in which the vorticity is not localized on lines, the analogy still works in detail — only, 
one should identify the current density with the vorticity field u. Table [T] summarizes the 
dictionary between magnetostatics and our system. 

In conclusion, the interaction between vortex lines in an incompressible fluid can be 
described exactly as magnetic interaction between currents. The only novelty is how the 
lines respond to these magnetic-type fields: given the peculiar structure of the lines' kinetic 
action in (94) (the first line), there is no analog of the Lorentz force. Instead, as discussed 
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Magnetostatics 


Incompressible Hydro 


— * 

current J 
magnetic field B 
vector potential A 


vorticity uj 
velocity field v 
hydroplioton A 



Table 1: Schematic dictionary between magnetostatics and our system of vortices in an incom- 
pressible fluid. 



in the previous section, the instantaneous velocity of a line-element is the local value of the 
'magnetic' field. 

We can now assess very easily the importance of our sound- mediated potential (66), 
relative to the purely kinematical long-distance interactions which, in our new language, are 
mediated by A. For instance, two vortex rings of radii -Ri,2, circulations ri^2, and orientations 
ni 2, at a large distance from each other have the same interaction potential energy as two 
magnetic dipoles, 

Vdip = ^ [3(f ■ /ii)(f ■ /i2) - /ii ■ /I2] , (97) 



with dipole moments 
We thus get 



fin = vr(r„i?2) f,^ . (98) 



Kiip ~ ^ {r,Rl) {T2RI) , (99) 



which is indeed a factor of (c^/f)^ bigger that the sound-mediated potential, eq. (75). Notice 
however that for very thin vortex cores, like e.g. in superfiuid vortex lines, eq (75) gets 
enhanced by a factor of (log i?/a)^, which can partially compensate the {v/csY suppression 
factor. 

— * 

The reader might be skeptical about the usefulness of our introducing the local field A to 
describe vortex interactions. In fact, aren't we supposed eventually to solve all the equations 
of motion? If we first solve the equation of motion for A and we plug the solution into the 
others, we are effectively reproducing the equations of motion deriving from (83). So, why 
bother introducing A in the first place? Although this is technically a valid viewpoint — no 
information is added by introducing A — such a statement is technically and conceptually 
equivalent to claiming that it is useless to introduce the local fields E and B (or V and 
A) in electrostatics and magnetostatics, since one can do everything at the level of the 
Coulomb force between charges and the magnetic force between currents. A position that, 
in hindsight, few would defend. We believe that, like for electrostatics and magnetostatics, 
our local rewriting of vortex interactions will prove valuable in the study of vortex line 
systems. 



7.1 Example: Kelvin waves 

To prove the usefulness of our approach, we now study via standard field-theoretical tech- 
niques the low-frequency spectrum of small perturbations of an infinite straight vortex-line. 
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Let's align the (unperturbed) line with the 2;-axis. The perturbed fields are 



X(A,t) =Xo + 7f(A,t) 
A{x,t) = Ao + T^{x,t) 



Xn = zX 



Ao = -z—\og{r/a) 
Zn 



(100) 
(101) 



where Xq and Ao are the unperturbed configurations, r = -^/x^ + j/^ is the distance from 
the ^-axis, and a is a UV cutoff (the line thickness). We can now expand Siocai in powers 
of vf and ijj. The linear term in the fiuctuations is, of course, going to vanish as Xq and 
Ao are solutions to the equations of motion. We are thus left with, at lowest order in the 
fiuctuations, a quadratic action: 



S 



dtdX \TtijTi\dt-K{ - j d^xdt \ {d'lp^y (102) 
+r j dtdX (^PiO) ■ dx^ + V7/;"(0) ■ vf + V4(0) • 7f Satt'I + ^^^^jAl{0) 



where '_L' and '||' are defined w.r.t. the 2-axis, and we have freely integrated by parts and 
thrown out boundary terms. In the second line, the notation ip{0) and Ao{0) reminds us 
that these 'bulk fields' and their derivatives are to be evaluated on the unperturbed line, 
at X = y = (the last three terms come from expanding their arguments in powers of if). 
One can thus see that the second-to-last term vanishes as O^Aq is odd in either x oi y (or 
simply vanishes for j = 3). Continuing further the decomposition of ip and vf into parallel 
and perpendicular components, and noticing that Ao only has derivatives in the orthogonal 
direction, we finally get 



dtdX \Teij-K\dt-K{ - / d^xdt \ {d'ip{y + T I dtdXip±{0) ■ dx^± 



S ~ —Wo 

- I d^xdt\ (5V||)' + r / dtdx (v±^/'||(o) ■ + \didfj^l,{Q) 



(103) 



Recall that we have not fixed the gauge yet. We would like to choose X = z = X^, that is, 

7r||(A,t) = 0. (104) 

Since vrH has completely disappeared from the action, we can perform this A — )■ ^ gauge choice 
directly at the level of the action without worrying about losing tt'I's equation of motion. 
Varying (103) with respect to tp^, ■i/'ii, and 7r_L we generate the three equations of motion: 









vV± + r<52(5'±)5.vfx 
vVii -rv±(52(f±) -Tfx 

dtfii. xz- 9,^x(0) + Vx^ii(O) + V±d^AliO) n] 



(105) 
(106) 
(107) 



These equations can be recast into a slightly more convenient form by insisting on manifest 
translational invariance: The background configurations Xo(A,t), Ao{x,t) (spontaneously) 
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break translations in the plane orthogonal to our line. As a results, these translations are 
non linearly linearized on the associated perturbations: 



7f_L ^ TT^ + el , V|| ^11 - ^± ■ V±aJ , (108) 

implying in particular — along the lines of Goldstone theorem — the absence of a gap at zero 
momentum. Although the action and equations of motion above are invariant under these 
transformations, the absence of a zero-momentum gap is not obvious at all, and will require 
non-trivial cancellations, since there are terms in which 7r_|_ appears without derivatives. It 
would be nicer to have a reformulation of these equations (and of the action), in which the 
broken translations act simply as a shift on vf and nothing else. This is easily accomplished 
by redefining the -^n field as 

i^\\=i^-^±- , (109) 

so that, using V5_^q = —T5'^{x±), the eom become 

= VV± + r52(xx)a,7f^ (110) 
= V^^ -dl^^-V^J^l (111) 
= dt^^ X z- d^^p±{0) + V±ip{0) . (112) 

Now 7rj_ always appears with derivatives, and translational invariance is realized simply as 

TT± 7r_L +^1. 

These equations describe waves propagating in the z direction. To see this, let's go to 
Fourier space for z and t (but not for x±), upon which the equations of motion become 

= {Vl-P)^± + r6\x^)ik7r^_ (114) 
= (yl-k^)^ + k^jf±-V±Al (115) 
= -iun^x z-iktf±{0) + V±i!{0) , (116) 

where now 7i± depends on k and u, and iIj± and ip depends on them as well as on x±. 

The shape of the perturbed line is uniquely determined by vf^. Therefore, as long as we 
are interested in just that, we only need the solutions for ipx and ip close to the line — see 
the third equation. Treating then V^^ in the first two equations as being of order 1/r^ ^ k"^, 
close to the line we get 

- r 

Vi^± + r52(fj_)ifc7f_L ~ ik\og{kr)7f± (117) 

27r 



Vi^ k^ vfx ■ ^ ~ ~ log(fcr) rl ■ vf^ (118) 

271 Att 



^The same change of variables at the level of the action, replaces the second line in ( 103 ) with 



J d''x[^{d,i>y -d,i>d,7T^-W^Al + ^{d,7r^-W^AlY] , (113) 



which has manifest shift-invariance for t:±. Varying the new action with respect to the fields one gets, after 
stratighforward manipulations, precisely the new equations of motion. 
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(The exact solutions at all distances involve a Bessel function of the second kind, KQ{kr), 
and have the precisely these asympthotics for kr ^ 1.) Plugging these solutions into the 
third equation, we finally get a wave equation for Tf± only: 

r 

- iu (vf^ X z) + — log(l/A;a) vf ^ = , (119) 

where we have stopped the limit r — t- at some UV cutoff a, on which, like before, the log 
weakly depends, and we have ignored a piece proportional to r]_d-\_ log(/cr)|^, which is finite 
for a — )■ 0, and can thus be thought of as redefining the UV cutoff a inside the log. Notice 
that we have expressed the log so that it is positive, since /ca -C 1. 

Such a wave equation describes two transverse, circularly polarized modes of oscillation, 
with dispersion relations 

CO = ±—P\og(l/ ka) , (120) 

and polarizations ;^(l,=Fi). Note that both these modes circulate (at a fixed height along 
the vortex line) in the opposite direction of the flow circulation. These are the famous 
Kelvin waves, and our results match their qualitative and quantitative properties as derived 
via Victorian methods [28] . 

Despite the (weak) dependence on the UV cutoff, the dispersion relation above has pre- 
dictive content. We can rewrite it in an 'RG-invariant' fashion by comparing the frequencies 
associated with two different momenta: 

or equivalently 

d(uj/k^) r , , 

— - = . (122) 

dhgk 477 ^ ^ 

Such relations are completely independent of the UV details of the vortex line, and are thus 
a robust prediction of the low-energy effective field theory. For instance, one can use them 
to measure indirectly the circulation F. 

7.2 Reintroducing sound 

It is now time to go back to the first part of our paper, and look again at how compressional 



modes couple to incompressional fluid flows. Consider first the lowest-order coupling (29) 



and, to be consistent with this section and the previous one, let's neglect the second term. 



which is a relativistic correction. In place of Vq, we could plug in the integral expression (79), 
which gives the incompressional velocity field as a function of the vortex lines' geometry. 
Clearly, the resulting Lagrangian term would correctly describe the interaction between 
vortex lines and sound, and can obviously be extended to all orders in (f/c^)^, following the 
expansion of Appendix |Xj The downside of this approach is once again the non-local nature 
of the resulting action: for each power of Vq, we end up with a line integral over all vortex 
lines. 
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Figure 4: Lowest order Feynman diagram contributing to the emission of sound by the scattering 
of two vortex rings. The dashed lines represent the hydrophoton and the wavy line the sound. 



The hydrophoton provides us with a conceptually much cleaner solution. Following the 
magnetostatics analogy, which, for the incompressional part of the fluid flow, is exact, we 
simply have 

Vo = V X A. (123) 

So, the local action that describes vortex lines, the hydrophoton, sound, and their interac- 
tions, is simply 



S =wo 



(124) 



+ / d^xdt {diAjY + \cl{V ■ 



+ ^ r„ /" dtdX dxXn ■ A{Xn, t)+ I d^xdt (V x Ay{{V x A) ■ V)ip' + . . 

The first line contains the vortex lines' peculiar kinetic terms. The second line is the free 
action for the hydrophoton and sound fields. The third line contains all the interactions: 
of A with the vortex lines, and of A with sound. The dots denote all higher order sound 
interactions of Appendix LaJ which, given the replacement (123), involve more powers of if) 
and A, and are straightforward to write down explicitly. 

The message is clear: vortex lines interact directly with A; sound interacts directly with 
A] vortex lines and sound interact only indirectly, via A. For processes that are amenable to a 
perturbative treatment, one can then use standard field theoretical/diagrammatic techniques 
to compute observables. Consider for instance a soft scattering process in which two vortex 
rings pass by each other at large impact parameter (compared to their radii). What is the 
sound emitted by such a "collision"? The most relevant Feynman diagram is depicted in 
fig. |4j the two vortex rings exchange a long-distance hydrophoton, and this internal line 
emits a phonon, as dictated by the AAip vertex in the action above. 



8 Discussion and Outlook 

We have initiated a systematic application of effective field theory techniques to the study 
of near incompressible hydrodynamical systems. As usual for local field theories, questions 
that can be dealt with in perturbation theory are reduced, at least conceptually, to "turning 
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the crank" : one can organize each observable as an expansion in Feynman diagrams — of ever 
increasing complexity as one moves to higher orders in perturbation theory — each of which 
can be computed in standard ways starting from the action. 

The small expansion parameters in our analysis have been the typical fluid-flow speed, 
which we took by assumption to be much smaller than that of sound, and the typical size 
of the regions where vorticose flow is localized, which we took to be small relative to the 
sound waves' wavelengths and to the distance to other vorticose regions. The smallness of 
the former parameter allows for a perturbative expansion in compressional modes (sound). 
The smallness of the latter allows for a multipole expansion. 

Our insisting on the field theory language has also led us to an entirely local action for 
vortex-line dynamics — which is nontrivial, given that, as is well known, each line carries a 
long-distance 1/r velocity profile with it. We have 'integrated in' such a velocity profile — or 
better, the associated vector potential (which we have dubbed the "hydrophoton" ) — thus 
ending up with a local field theory describing vortex lines, the hydrophoton field, sound 
modes, and their mutual interactions, to all orders in {v/cs). In this reformulation, vortex 
lines interact with themselves, with one another, and with sound, by exchanging hydropho- 
tons. We believe this to be substantial technical improvement over the standard "vortex 
filament" model. Once again, questions amenable to a perturbative expansion are straight- 
forward to analyze — and answer. 

From a theoretical viewpoint, there are a number of open conceptual questions about our 
field-theory formulation of vortex line dynamics. First, we did not construct it bottom-up 
from standard low-energy effective field theory principles, nor did we derive it top-down 
from our general fluid action ^ — we just guessed the vortex-line action, and checked that it 
reproduces the correct dynamics as derived from the Euler equation. To be responsible field 
theorists, we should remedy this. In particular, as we emphasized, the dynamics of vortex 
lines are strikingly different from those of string-like objects in empty space. Vortex lines 
obey first order equations of motion, which means that their local velocity is completely 
determined by the geometry of the lines' configuration. According to standard effective field 
theory intuition, all robust properties of the dynamics should follow from symmetry — in order 
not to be disrupted by renormalization, for instance. So, what is the symmetry protecting 
such a peculiar behavior? Compared to string-like objects in empty space, we seem to have 
less symmetry — Lorentz or Galileo boosts are gone, being (spontaneously) broken by the 
surrounding medium. This allows us to write down more Lagrangian terms, including our 
peculiar kinetic term with one time derivative — but the action (83), or equivalently (94), 
is hardly the most generic action invariant under translations and rotations. This is all 
the more frustrating as our original fluid action, eq. ([6]), was constructed precisely as the 
most generic action invariant under all the appropriate symmetries. Reparameterization 
invariance under A — X{X,t) — which is symmetry of (83) and (94) — will play of course a 
crucial role in reducing the number of allowed Lagrangian terms. We plan to address these 
symmetry questions in the near future. 

Second, for small vortex rings we would like to go one step further, and treat them as 
point-like objects with multipole moments (starting with a 'magnetic' dipole), and a number 
of degrees of freedom: their position, radius, orientation, and an infinite tower of excited 
states (Kelvin waves). Once again, given the first-order nature of the equations of motion, 
the dynamics of these objects are quite peculiar. For instance, from a preliminary analysis 
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it seems to us that the orientation is not a low-energy degree of freedom: one can certainly 
choose any orientation as an initial condition, but then one cannot change it via low-frequency 
processes. In a sense, unlike for ordinary string loops in empty space, the orientation degree 
of freedom seems to belong in the excited Kelvin wave spectrum. Again, all the vortex 
ring properties should follow from systematically applying symmetry considerations to the 
effective theory of point-like objects in the fluid, along the lines of [8j. 

Third, what is this hydrophoton anyway? We have a local field A, whose equation of 
motion in the incompressible limit, 

\/'^A = u, (125) 

implies instantaneous signal propagation at arbitrarily large distances. This is clearly an 
artifact of the incompressible approximation, and once the speed of sound is brought down 
to finite values, we expect A to propagate at that speed — no small disturbances can propagate 
faster. This sounds very much like sound. Is A in some sense sound itself? 

But, while we responsibly address all the conceptual points raised above, nothing pre- 
vents us from looking for relevant physical applications of our methods, which can range 
for instance from classical vortex ring systems in ordinary fluids like water [21], to quan- 
tized vortex line interactions in laboratory superfiuids or in pulsars, where dense arrays of 
vortex lines in the neutron superfiuid are expected to interact strongly with much denser 
arrays of magnetic flux tubes in the proton superconductor, with crucial consequences for 
the dynamics of the star whole f3T]. 
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A Systematic expansion of the Lagrangian 

Here we make more explicit the steps taken to arrive at the foundation of the work in this 
paper. We have attempted to make this section as complete as possible, and as a result there 
is some overlap with the more rapid expansion found in Section [2| and with refs. ^T3[ [1]. 

We begin with the field theoretical description for a relativistic perfect fluid [H [131 112] • 
Its action is given in eq. (|6|, which we reproduce here for completeness: 

S = jdt^xF{h), with h = Vdet B^^ , (126) 

where B^^ = d^cp^d'^cp'^ and where F is a generic function. A particular functional form of F 
corresponds to a particular equation of state. However, in order carry out the expansion of 
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T\-l 



the Lagrangian in fluctuations about a background flow, there is a slightly more convenient 
starting point than eq. (126). Manipulating b we can pull out the Jacobian determinant for 
the equal-time x -H- mapping, which, for nearly incompressible flows, is close to one 

= det9^0^9^0^ 

= det(90^-90-0®(^ 

= det (d(f)'^ ■ (1 - {{d(j)'^y^ ■ 0) ® {{d(f) 

= {det dcj))"^ det {l - V ^ v) 

where we defined the matrix {d(f))ij = di(j)j, and the vector 

V = -{d(f)'^)-^ ■ $ . (128) 

Notice that, by the implicit function theorem, v is precisely the usual fluid velocity field — 
hence the name: 

dx{(f),t) 



(127) 



dt 



(129) 



The second determinant in (|127|) is easy to compute by going to a basis where, locally, 

(130) 



the X axis is aligned with v. We get 

b^ = (det90)2(l- I^TH 
and therefore we can rewrite the Lagrangian in the form: 

C = -woc^f ( det d(f) - \v\^ ) 



(131) 

where / and are defined in eqs. ([12]), ([13]). As <\){x,t) is one-to-one with x for all t, we 
can change coordinates to those of the comoving volume elements and we have arrived at 
equation ([O]): 



S 



j d^(f)dt det J /((det J^i) v^r^Wc^) • 



(132) 



J*j is the Jacobian matrix dx'^/d(f)^, and v = dtx{(f),t). 

We now expand this action in powers of v/cs — which is small by assumption — and of the 
compressional field if, which, as explained in sect, [ij describes the small deviations from a 
purely volume-preserving time evolution: 



f(0,t) = fo(0,t) +ljj{(l),t) 



det ( Jo) = det 



dxQ 



(133) 



For V and det J the expansion in ip truncates at finite order (the determinant of a 3 x 3 
matrix is a cubic function of that matrix): 



v = vo + dtip 



vo = dtXQ 



(134) 



det J = 1 + [Vo^] + ^ ([Vo^]' - [Vo^']) + I {[Vo^f - 3[Vo^][Vo^2] + 2[Vo^']) , (135) 
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This separation is also convenient for expanding around a static background, (f>^ — , and was utilized 



in p. 
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where [...] means the trace, VqV'" means the n-th power of the matrix dip'^ /Bxq, and we have 
used the identity 



Of course, for the inverse determinant and the square root entering the action (132) the 



expansion does not truncate, and goes on to all orders. Up to cubic order in ip we get 



S = wo I d'<f)dt \+lv' + i^^^ (137) 



^ (-3c,^[Vo^] ((5,^)^ - c2[Vo^']) + cVa [Vo^]^) + • • • } 



We have used /'(I) = 1, /"(I) = c^/c^, /s = /'"(I), and for notational simplicity we 
have dropped the subscript zero in vq, and so, from now on, f = fo is the underlying 
incompressional velocity field. Even though /s is a free parameter, as pointed out in [1] one 
typically expects 

fs - cl/c' . (138) 

We have dropped total derivative terms, most notably the term —wqc^ /(I) det J, since det J 
itself is a total derivative — schematically: 

det J = det d^x oc ee ■ dff,x d(j,x d^x = d^j, ■ [ee ■ x d^x d^x) . (139) 

We have also dropped terms that are total derivatives with respect to xq, even though the 
integral is in d^(f), for reasons that will become clear in a moment. Notice that for different 
powers of ^Z', we have stopped the expansion at different orders in vq. The reason is that the 
two expansion parameters are not entirely independent, as we will explain below. Notice 
also that the first line is just the action for the incompressible fluid flow, in the absence 
of compressional perturbations. To be used in this sense, it should be supplemented by a 
volume-preserving constraint for xo{4>,t), as explained in the main text. However that is not 
our goal, since we are interested in solving for the dynamics of ip in the presence of a given 
background incompressible fluid flow Vq. We will therefore discard the first line for what 
follows. 

We are not finished yet. We want to now go into the xo coordinate system. The reason is 
that in an experiment it is much more convenient to parameterize a velocity field in terms of 
its dependence on the physical x coordinates rather than on the comoving ones. Since the x 
and Xq coordinate systems coincide at lowest order in ip, using xq will suffice for lowest order 
computations. This is an easy transformation for most things in the action: The Jacobian 
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Jo has unit determinant, and all the spacial partial derivatives are with respect to xq already 

(this is the reason why above we were allowed to neglect total ,Xo-derivatives). The non- 
trivial pieces are the partial time derivatives as they are taken at constant (j) values. We can 
write the time derivatives in a more compatible form: 

dt^ixoicf), t),t)\f = dt^ixo, OUo + vi^rixo, t) . (140) 

OXq 

— * 

So, changing coordinates from to and inserting the expanded time derivatives in the 
action above we arrive at an action which can be organized in the following way: 

where the n's above mean 'positive powers of f' (1, 2, 3, etc.). Dropping for notational 
convenience all the subscript O's, and denoting partial time-derivatives by overdots, we have 
explicitly: 

V ^woj d'xdt ^ - clidiP]'') (142) 

V -n^oj d'xdt ["^mmf] + "-^m' - ^j'm} (143) 



.n^woj d^xdt ■ if+Viiv ■ - ^ v^[djlj] 
+ ^""'f^ {f^ ■ ^+ ■ ^W) + ^ AQi>\ + • • • } (144) 

d^xdt {0,(1; • v)i/y' - '^{v-i))\d^\ + 2 ((^ ■ ^Wf 



c 



ij (^ - c9 (v ■iif + ^(c'- cj) + ■ ■ ■ } (145) 



While all of the above looks like a total mess (and of course there is an infinite tower of 
terms), it is both easy to generate and easy to interpret perturbatively. In particular: 

• S^2 describes the free propagation of compressional modes (sound waves); Cg is indeed 
their propagation speed. The associated (Feynman) propagator is 

(T(V^V)) = 1— ^ . (146) 

Wq UJ^ — Cjp ''+16 

• S^s describes the sound waves' trilinear self-interactions; in a Feynman diagram these 
will correspond to a vertex with three lines attached. 
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• yn corresponds to 'tadpole' diagrams for ip — diagrams with a single ip line attached 
to the external source vq: In the presence of a non-trivial background incompressional 
velocity field, = is not a consistent solution. Terms in the action that are linear 
in ip describe, at lowest order in ip, how a non-trivial ip field is generated by such a 
background velocity field. 

• 5*^2 corresponds to interaction vertices where two ip lines attach to the external 
source vq, describing for instance the scattering of sound waves by vorticose motion, 
as in sect. HI 

• And so on. 

As discussed throughout the text, note that the terms with the additional c^/c^ factors scale 
like many of the others in terms of powers of v/cg, derivatives, and so forth. These terms 
are the relativistic corrections, and will be necessary to correctly describe ultra-relativistic 
fluids. 

By comparing 5*^2 and S^^r, we see that, for small v/cg, the ip field generated by the 
background fluid flow scales like ip v'^ (we are implicitly assuming that v scales like f^, 
keeping the typical length scales fixed.) According to this power-counting scheme, the terms 



explicitly displayed in (142)-(145) are all the terms in the action up to order v . Things 
become more complicated when not all of the ip^s appearing in a diagram are "generated" 
by the background Vq, like for instance in a scattering process. Also, in general there are 
important power-counting differences between highly off-shell internal lines with u <^ Cgp, 
and on-shell external lines with u = Cgp. These subtleties for power counting in classical 
perturbation theory are addressed systematically for general relativity in ref. [8J. It would 
be interesting to derive analogous power counting rules for our case. We leave this for future 
work, since for our lowest-order computations in this paper a quick case-by-case analysis 
suffices. 



B Scattering cross section and emission rate from the 
amplitude 

For the benefit of the reader we derive the scattering cross section formula for one-to-one 
scattering. Following Srednicki's [19] quick construction — which avoids the use of wave 
packets — we start with the probability of scattering: 







2 




l/> (' 









> 


(/|/>(' 


;|0 



p = I I = . (147) 



Note that our M here is different then the usual one because we have not yet removed any 
delta-functions. Assume that we are performing our experiment in a large box of volume 
V and over a time T. The norm of a single particle state is given by (with the so-called 
relativistic normalization, which is in fact still convenient for non-relativistic systems [1]) 

{k\k) = 2E {27V f 6^0) =2EV . (148) 
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For our 1 — )■ 1 scattering off external sources of sect.|4| let's take the incoming momentum 
and energy to be pi and Ei, and the outgoing ones to be p2 and E2. Additionally, we want 
to sum over final momentum states We then have 

which is in the continuous limit becomes 



To get the total cross section we need to divide by the incoming flux (which for an 
incoming particle moving at speed Cg is simply [cg/V)) and divide by the total time T. 
Finally, we arrive at the general expression 

1 1 d^p2 \Mipi,P2)\^ 

Ai can be computed following the standard relativistic Feynman rules [1], with the caveat 
that whenever an external source J{x) appears in a Lagrangian term, its (four-dimensional) 
Fourier transform should appear as a factor in A^. As a consistency check, notice that 
for Lagrangian terms with no external sources, we formally have J{x) = 1, whose Fourier 
transform is the usual (27r)'' (5(energy)5^(momentum). 

Now, if our vorticose source is time-independent, then we have an energy-conserving delta 
function in our A^, such that 

M M'{27v)6{Ei- E2) (152) 
\M\^ \M'\\27i)6{Ei- E2){27t)S{0) (153) 
= \M'\\27v)6{Ei- E2)T (154) 

and thus, in this time-independent limit we have that 

1 1 d^P2 I , ,|2, 

7g2E~i {2n)^2E2 



c^cr(i^i)(time independent source) = — — -7^r^^\M'iPi^P2)\\^^)HEi-E2) (155) 



which, as a good check, matches Peskin and Schroeder's Rutherford scattering problem in 
their fourth chapter |20j . 

Similar considerations apply to the emission process discussed in sect. |3} and yield 



eq. (31) 



C Scattering sound waves off sources: Checks 
C.l Discrepancy with previous results 

The results found in [TH] were developed in the context of classical scattering theory with 
the same set of assumptions as ours: in particular, low Mach number of the source flow, low 
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intensity and high frequency relative to the typical source frequency of the incoming sound 
waves. Moreover, the computations of ref. flS] are appropriate for non-relativistic fluids only, 
so for the purposes of this comparison we will set c^/c — )■ in our results. 

The first step in the comparison is to note that the vorticity field, defined a.s w = V x v, 
satisfies 



(pi X P2) ■ wiApf") = i (pi ■ Ap){p2 ■ v) - {p2 ■ Ap){pi ■ v) 



(156) 



where w{Ap'^), v{Ap^) are four-dimensional Fourier transforms, and for brevity we are using 
a relativistic notation, Ap^ = {Au},Ap)). Defining the scattering angle through the usual 
Pi ■ P2 = cos 9, we can write 



iPi X P2) ■ w{Ap^') =i{l- cos 9) pi [{pi +P2) ■v] + 



/ Acq 

\CsPl 



(157) 



Here Aco/csPi is small — see sect, [i] — and so we will neglect it in the following. 
Now, using (157) we can re-express (51) as 



da 



dnd{Auj) Act (2vr)3 T (l-cos^)^ 



cos^ 9 „ , _ , , ,,m2 / Auj 



CsPl 



(158) 



Using our notation for the Fourier transform (which is different than that of [IB]) we can 
express 



(159) 
(160) 



Wk{Ap^')w;{Apf') = I d^xd^y e-''^P'''^+'^P'y^wi{y^')wk{x^') 

d^R dA e<^-^-^>^) 



X / d^x 



dr wi{R + x,T + A/2)wk{x, r - A/2) 



where we have shifted the integral variables appropriately. Using the notation of [18] we can 
write our final scattering cross section as 



da (27r)z/2 cos^^ 



AkAi Ski{q, V - z/q) 



d'a2dv 4c4 (l-cose)2 

where v = 002, vq = uji, q = Ap, dVl2 = dVt, Aj, = {pi x p2)k 

1 



'IQV. 



Skiiq, y - vq) 



and where 



^ , rfV rfA eH^^-^-^f^J / rf3x(«;fc(i? + x,A/2)u;z(x,-A/2) 



(2vr) 



(/(A/2)/(-A/2)) ^ i drfir + A/2)/(r - A/2) . 



T 



(162) 
(163) 



T/2 



When we compare (161) to that of equation (19) in we see that there is a factor of 2 
discrepancy. 
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C.2 Agreement with optical theorem 

The optical theorem as expressed in quantum field theory is a consequence of the unitarity 
of the S'-matrix. In our case, it relates the imaginary part of the one-to-one amplitude in the 
forward scattering limit, to the sum over the squares of the one-to-anything amplitude^^ 



2 Im M{i -^i) = J2 J dUf \Mit f)f 



(164) 



It is a statement about full amplitudes, but perturbatively it can be expressed as a relation 
between particular graphs. To lowest order (in the coupling) the relationship is simply given 
by: 



2 Im 




{2tt)^2uj 




(165) 



We can rewrite the right hand side of the above equation as the one-to-one scattering cross- 
section using (151). Written this way, we have 

2 Im M{pi,pi) = 2uiCsT a^i^i-) . 



(166) 



Eq. (51) provides us with the right hand side of this equation. We need to calculate the left 
hand side. 

The imaginary part of the amplitude comes from the ie prescription in the internal prop- 
agator of sound modes. Using generating functional methods to ensure that all symmetry 
factors are correct, we have 



{0\T^\xi)tlj^{x2)\0)^ 



1 



6 



I 6J^{xi) i 5Ji{x2) 2! 



1 6 



X ■ 



3! 



d^xd^x' j\x)A''\x - x')J\x') 



Ml 



(167) 



Upon taking all the functional derivatives it is easy to see that one is left with four terms. 
Inserting this correlation function into the LSZ formula the amplitude is 



dujq d^q 



1 



yiy2 



2n {2ny^ — + c^g ^ — ie 

X [wi (%i) ■Q)+^q (%i) ■ Pi)] {q ■ Pi 
X [^2 {v{y2) ■q)+^q {v{y2) ■ P2)] {q ■ P2 



g«2/i(pi-i3)g«y2(g-P2) 



(168) 



where we have suppressed relativistic corrections for algebraic simplicity; their inclusion is 
straightforward. 



^We use the conventions given in ^IQj and [20], which should be noted are different than those given in 



36 



Taking the forward scattering limit, p2 Pi, and performing the y integrals which yield 
Fourier transforms of the velocity fields we arrive at 



iM{pi ^pi) 



dojq d^q 1 



(169) 



where = {q ■ pi) [wig* + ujqP\\- In this form it is easy to see that the only contribution to 
the Im Ai comes from the ie, we have 



Im M.{pi,pi) 



duq d^q TT 



271 (27r)3 2c,|gl " 
X v\pi - q)V{q - pijVV^ 



[5{ujq - Cs\q\) + 5{uq + Cs\q 



where we have isolated the imaginary part with the standard formula 

— -— = V- =F moix) , 

X ±lt X 



;i7o) 



(171) 



where V is the principal value. 

Only the first delta-function in (170) is going to overlap with the support of v^{j)i — q). 
This is because v^{pi—q) offers a narrow width around zero frequency relative to the incoming 
frequency (as the time scale of the velocity flow is much longer than the incoming frequency). 
For instance, in the static limit v'^{pi — q) oc S{ui — Ug)v\pi — q). For the same reason, as 
discussed in Section |4| to lowest order in the energy transfer we can take ui = U2 = uj. 
Evaluating this final delta function, relabeling our variables of integration and performing 
some trivial algebra we have: 



2lm M{pi,pi) = {2uCsT) 



dnd{Au) 
(27r)3 4^ 



(P1-P2)' e{Auj,Ap)-{p\+pi,) (172) 



Inserting this into the left-hand side of ( 166 ) and comparing with our scatting cross section 
(51) calculated in the bulk of the paper we see that indeed the optical theorem is satisfied. 
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